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Abstract 

The interaction between vacancies and edge dislocations in face centered cubic met- 
als (Al, Au, Cu, Ni) is studied at different length scales. Using empirical potentials 
and static relaxation, atomic simulations give us a precise description of this in- 
teraction, mostly in the case when the separation distance between both defects is 
small. At larger distances, elasticity theory can be used to predict this interaction. 
From the comparison between both approaches we obtain the minimal separation 
distance where elasticity applies and we estimate the degree of refinement required 
in the calculation. In this purpose, isotropic and anisotropic elasticity is used as- 
suming a perfect or a dissociated edge dislocation and considering the size effect as 
well as the inhomogeneity interaction. 
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1 Introduction 



Understanding the interaction between dislocations and vacancies is a key 
point in material science. In metal plasticity, this interaction controls disloca- 
tion climb and thus the plastic strain of a metal at high temperature when the 
glide of dislocations is blocked by obstacles like precipitates or the dislocation 
forest. This interaction also plays an important role in phase transformations. 
Indeed, as dislocations act as sinks and sources for vacancies, they control the 
way vacancies adjust their concentration in non-isothermal kinetics and thus 
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the diffusion. Moreover, the vacancy super-saturation along dislocations leads 
to an enhancement of diffusion in their vicinity: dislocations behave thus as 
short-circuits, a phenomenon known as pipe diffusion. 

The interaction between vacancies and dislocations can be modeled using elas- 
ticity theory (for a review, see Ref. 1). In 1949, Cottrell and Bilby calculated 
the so-called size contribution, showing that the interaction energy is equal to 
the pressure created by the dislocation times the relaxation volume due to the 
vacancy [2] . This energy is thus decreasing with the inverse of the distance be- 
tween both defects. In 1961, BuUough and Newman applied results of Eshelby 
theory [3,4] and brought out another contribution, the inhomogeneity interac- 
tion [5]. This contribution is proportional to the square of the pressure and of 
the shear stress created by the dislocation and is always attractive. It is much 
shorter range than the size contribution as it is decreasing with the square of 
the inverse of the distance. Such results are based on a continuum description 
of matter and assumes that a vacancy can be modeled by an inclusion having 
different elastic constants from the matrix. 

First atomic calculations of the interaction between vacancies and disloca- 
tions [6-11] were used to check the vahdity of the elasticity theory and to 
obtain results in the vicinity of the dislocation core where a continuum de- 
scription breaks down. For face centered cubic (fee) metals [8,11], the vacancy 
binding energy was found to reach a maximum in the region of the stacking 
fault rather than in the immediate vicinity of the partial dislocations. More 
surprisingly, the vacancy was found to be attracted to a dissociated edge dis- 
location both above and below its glide plane. Some authors interpreted that 
as a manifestation of the inhomogeneity interaction [8]. In the same way, for 
body centered cubic metals, an anomalous positive binding energy was found 
in some cases on the tensile side of the dislocation elastic field [6, 7] and the 
comparison with elastic models show discrepancies even for large separation 
distances [6,10]. Most of these simulations have been performed before the de- 
velopment in the mid-eighties of reliable empirical potentials for metals, like 
the embedded atom method (EAM) [12] or the tight-binding second-moment 
approximation (TB-SMA) [13, 14]. They were based on pair potentials which 
fail to correctly predict elastic behavior. In particular, all cubic crystals mod- 
eled with pair potentials have to obey the Cauchy relation, which restricts 
the elastic constants C12 and C44 to be equal. Therefore, one can legitimately 
wonder whether results obtained with these potentials are not biased. On the 
other hand, EAM and TB-SMA potentials do not suffer such restrictions and 
have been shown to be well suited to study dislocation behavior. Some au- 
thors used these potentials to study the behavior of vacancies in the vicinity 
of dislocations, in particular to simulate pipe-diffusion in fee metals [15-20]. 
These studies showed that the most stable position of the vacancy lies at the 
edge of the supplementary half plane corresponding to a partial dislocation 
and that diffusion is faster along dislocations than in the bulk. Nevertheless, 
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no careful comparison between atomic simulations and elasticity theory was 
performed. Such a comparison is useful to determine the range of validity of 
a continuum representation and the degree of refinement needed by elasticity. 
Indeed, not all materials properties can be modeled at the atomic scale and 
a description of the behavior of the material intrinsic defects at a larger scale 
is still needed. It is particularly true for plasticity where properties are linked 
to the collective behavior of dislocations which has to be modeled at a larger 
scale than the atomic one. It is thus interesting to see what conclusions can 
be drawn at the atomic scale with empirical potentials and how these results 
compare with predictions of the elasticity theory. 

In the present work, we study the interaction between vacancies and dis- 
locations in fee metals at the atomic scale using EAM potentials and at a 
mesoscopic scale with the elasticity theory. We hmit the discussion to edge 
dislocations as creep is controlled by their climb and thus by their interaction 
with vacancies when an external stress is applied. Several metals having the 
fee structure are studied so as to bring out general results for fee metals and 
identify a possible influence of the material properties, like the stacking fault 
energy or the elastic behavior. The results presented in this article are directly 
obtained from atomic simulations and then compared with predictions of the 
elasticity theory using different levels of approximation. 



2 Atomic simulations 

2. 1 Interatomic potentials 

The atomic simulations in this work are performed using empirical potentials 
of the EAM type. Such potentials are widely used to study dislocation proper- 
ties as they present a good compromise between the size of the system which 
can be simulated and the realism of dislocation simulations. So as to test the 
generality of the results, various metals with the fee structure are examined. 
The corresponding potentials are the ones developed by Ercolessi and Adams 
for Al [21], by Medlin et al. for Au [22], by Mishin et al. for Cu [23], and 
by Angelo et al. for Ni [24]. Some potentials leading to a better description 
of these metals properties may exist in the literature, but the purpose of this 
study is mainly to draw the main features of the vacancy-dislocation inter- 
action which are common to fee metals and to test the validity of the elastic 
theory compared to atomic simulations. In this purpose, the chosen EAM po- 
tentials appear to be well-suited as they cover a wide range of properties, as 
summarized in Tab. 1. For instance, Au potential leads to a low stacking fault 
energy 7 and Al potential to a high one. Concerning their elastic behavior, 
Al potential appears to be almost isotropic ~ 1) whereas Cu potential is 



3 



strongly anisotropic. Moreover, Al potential gives quite soft elastic constants 
whereas Ni potential gives hard ones. At last, the vacancy energy of formation 
obtained in Al is low compared to the one in Ni. 

2.2 Simulation cell 

The main characteristics of the unit cell used for atomic simulations are 
adapted from those used in [26-28] and are sketched in Fig. la. The x, y, and 
z directions are parallel, respectively, to the [110], [112], and [111] directions. 
With such orientations, perfect edge dislocations can easily be introduced with 
a glide plane lying between two (111) planes, a line direction along the [112] 
axis and a Burgers vector along the [110] axis. As periodic boundary condi- 
tions are used in all directions, two edge dislocations with opposite Burgers 
vectors have to be introduced in the simulation box. The distance between the 
two corresponding glide planes is chosen as half the height of the simulation 
box. According to elasticity theory, these two interacting dislocations form a 

Table 1 

Properties of the different EAM potentials: lattice parameters a, vacancy energy 
of formation E\j, vacancy relaxation volume (5J7v) stacking fault energy 7, elastic 
constants Cn, C12, and C44, elastic anisotropy factor A, isotropic shear modulus ji 
and Poisson coefficient v given by Voigt average (Eq. 3), and dissociation distance as 
given by isotropic (o!'^°) and anisotropic (c?^°') clastic theories and atomic simulations 



Potential 


Al [21] 


Au [22] 


Cu [23] 


Ni [24,25] 


a (A) 


4.032 


4.079 


3.615 


3.52 


E{, (eV) 


0.693 


0.873 


1.272 


1.592 


<5Qv (A3) 


-6.6 


-10.0 


-3.5 


-1.9 


7 (mJ.m~^) 


104 


30.7 


44 


89 


Cii (GPa) 


118 


185.8 


170 


246 


C12 (GPa) 


62.3 


157.1 


123 


147 


C44 (GPa) 


36.7 


38.9 


76.2 


125 


A = 2C44/(Cii - C12) 


1.32 


2.71 


3.2 


2.5 


H (GPa) 


33.2 


29.1 


55.2 


94.7 




0.320 


0.418 


0.324 


0.277 


(A) 


11.7 


43.4 


37.0 


27.5 


d'^" (A) 


11.7 


41.7 


35.0 


26.7 


(A) 


15 


40 


31 


25 
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stable dipole for a 7r/4 angle between the vector linking the centers of the 
dislocation lines and the glide planes [29]. This result holds for an isolated 
dipole, but, due to periodic boundary conditions, the dipole in the simulation 
box interacts with an infinity of dipoles. We thus choose to put a distance 
equal to the half of the unit cell dimension in the [110] direction and we find 
this position to be stable. One should notice that dislocations in the simulation 
box are straight dislocations and do not have any kink nor jog. 

The atomic positions and the shape of the unit cell are relaxed at K using a 
conjugate gradient algorithm. As expected, each edge dislocations dissociates 
into two partial dislocations separated by a stacking fault [29,30]. Atoms be- 
longing to the core of these partial dislocations do not show fee nor hexagonal 
compact neighborhood. They can therefore be easily identified in the simula- 
tion box (Fig. lb), allowing the "measurement" of the dislocation separation 
distance d'^^". In agreement with the variations of the stacking fault energies 7, 
this separation distance is small for Al, big for Au and in-between for Cu and 
Ni (Tab. 1). The convergence of the dissociation distance d^^° with the dimen- 
sions of the unit cell has been checked and is found to be good for stackings 
of 160 and 60 planes respectively in [110] and [111] directions 



2.3 Vacancy- dislocation interaction 

Using the different EAM potentials, we calculate the vacancy-dislocation in- 
teraction energy for all possible positions of the vacancy in the simulation box 




Fig. 1. (a) Sketch of the simulation cell showing the two dissociated edge dislo- 
cations and their glide plane, (b) Aluminum simulation box containing two edge 
dislocations. For clarity, only atoms belonging to the dislocations are shown: atoms 
in yellow (clear) are in the core of the partial dislocations and atoms in red (dark) 
in the stacking fault. The cell contains respectively 160, 120 and 60 in x, y and z 
directions, corresponding to 190800 atoms. 
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Fig. 2. Vacancy-dislocation interaction energy Ey^^^ in (a) Al, (b) Au, (c) Cu, and 
(d) Ni as given by EAM potentials. On the left part of the figure, only attractive 
positions for the vacancy are shown (Ey^^ < 0), whereas on the right part only 
repulsive positions are shown (Ey^^ > 0). The positions of the partial dislocations 
and of their glide planes are sketched by lines. 
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Fig. 3. Vacancy-dislocation interaction energy Ey^^ in (a) Al, (b) Au, (c) Cu, and 
(d) Ni as a function of distance from the center of one of the two partial dislocations. 
Scans are performed in the [110] direction for different (111) planes lying above (in 
the compressive zone) or below the glide plane (in the tensile zone). The origin is 
defined as the center of one of the two partial dislocations and axes are defined on 
Fig. la. 



Table 2 

Lowest vacancy-dislocation interaction energy Ey^f^ (in eV) corresponding to a 
vacancy located on the edge of a partial dislocation. 



Authors 


Al Au 


Cu Ni 


Present work 


-0.181 -0.190 


-0.261 -0.279 


Huang et al. [15] 




-0.37 


Huang et al. [16] 




-0.27 


Hakkinen ct al. [17] 




-0.14 


von Boehm et al. [18] 


-0.11 


-0.22 



with respect to the dislocations. This energy is defined as 

-E^V-D = E + Ey-D — E-D — Ey, (1) 

where E, -Ev-Dj E^, and Ey are the energy of the simulation box without any 
defect, with both the vacancy and the dislocations, with the dislocations only 
and with the vacancy only, respectively. All these energies are calculated after 
relaxing the atomic positions and the shape of the unit cell. 

The calculations results are displayed both as functions of the vacancy po- 
sitions (Fig. 2) and as energy profiles along lines parallel to the glide plane 
(Fig. 3). As expected, compressive regions attract the vacancy whereas tensile 
zones are repulsive. For potentials presenting a low stacking fault energy, an 
attractive and a repulsive basin can be associated with each one of the partial 
dislocations, as they are well dissociated. For Al it is not so clear since the 
separation distance between the two partial dislocations is smaller. 

For all potentials, the most attractive position is located in the core of the par- 
tial dislocations, on the edge of the supplementary half-plane above the glide 
plane in the compressive region. This is in disagreement with first results ob- 
tained with pair potentials [8,11] but confirms most recent studies using most 
rehable empirical potentials like EAM [17], TB-SMA [18], or pseudopoten- 
tials [15,16]. The corresponding interaction energies are displayed in Tab. 2 
and compared with previous calculations, showing a dependency on the po- 
tential used for a given metal. It seems that the higher the vacancy formation 
energy, the stronger the binding to the dislocation. Using the same potential 
for Al, Hoagland et al. [19] calculated the interaction energy of a vacancy with 
a screw dislocation {E^^^ = —0.188 eV) and Picu and Zang [31] with a 60° 
dislocation {E^^'^^ = —0.163 eV). The values obtained in these works are close 
to the one we find for an edge dislocation (£^v*-d — —0.181 eV). Therefore, 
the binding of a vacancy to a dislocation seems to be quite independent on 
the dislocation character. 

Looking at the variation of the interaction energy in the plane just above the 
glide plane in the compressive region (Fig. 3), one can see that the binding in 
the stacking fault is smaller than the one in the core of the partial dislocations. 
We do not observe a repulsion in the middle of the stacking fault as observed 
by Hakkinen et al. [17]. Thus a Suzuki segregation of vacancies is expected. 
Nevertheless, this large difference of the binding energies indicates that va- 
cancies cannot easily migrate from one partial dislocation to the other as they 
have to cross the stacking fault. In such a case, one expects that diffusion is 
faster along the partial dislocation than in-between. This was shown by Fang 
and Wang [20] who used the same potential to calculate migration energies 
of a vacancy in Al. They obtained a higher energy for transverse core diffu- 
sion than for longitudinal core diffusion. Therefore, fast diffusion should occur 
along two well-separated pipes corresponding to the partial dislocations. Pipe 
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diffusion in a larger zone corresponding to tlie two pipes linked by a ribbon, 
as observed by Huang et al. [15, 16] in their molecular dynamic simulations, 
should only occur at high temperatures, where the partial dislocation cores 
spread over the stacking fault. 

Above the glide plane, in the tensile zone, most of the vacancy positions are 
repulsive and only a few ones are slightly attractive. This contrasts with first 
atomic calculations which were performed with pair potentials and predicted 
an attraction of the vacancy above and below the glide plane [8, 11]. 

Concerning the configuration in which the repulsive interaction is maximum, 
one cannot draw a general situation. Indeed, for An and Cu this position is 
located above the core of the partial dislocations, whereas for Al and Ni it is 
located in-between. Moreover, for Al and Cu, it hes three planes below the 
glide plane, for Au two planes below and for Ni only one plane below, i.e. 
in the habit plane of the stacking fault. We try to rationalize these differ- 
ences by looking at the interaction energy Ey^^p between the vacancy and 
the stacking fault (Fig. 4). When the vacancy lies in the habit plane of the 
stacking fault, this interaction is repulsive for Al and Ni, repulsive too for 
Au, but an order of magnitude lower, and attractive for Cu. This opposite 
interaction in the case of Cu arises from the fact that the stacking fault in 
Cu is in compression, whereas it is in tension in Al, Au and Ni. This was 
seen by computing the atomic pressure from the EAM potential using the 
definition of Vitek and Egami [32] (Fig. 4). This repulsive interaction with the 
stacking fault in Al and Ni adds, in the tensile zone, to the elastic repulsive 
interaction with the partial dislocations. Therefore the strongest repulsion is 
found between the partial dislocations. On the other hand, in the case of Cu, 
there is a compensation between the interactions with the stacking fault and 
the partial dislocations. Nevertheless, this explanation which tries to spfit the 
interaction in two different contributions does not really allow to predict the 
binding or the repulsion of the vacancy close to the dislocation. This shows 




-4 -2 2 4 

Distance in [111] direction (dm units) 




-4 -2 2 4 

Distance in [111] direction (d^ j | units) 



Fig. 4. Vacancy-stacking fault interaction energy -Eyl^gp (left) and atomic pressure 
P (right) as functions of distance from the stacking fault. Distances are normalized 
by dill, the interplanar distance in [111] direction in each material. 
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the usefulness of an atomic calculation close to the dislocation cores where the 
vacancy behavior is hard to predict. 



3 Elasticity 



As long as the vacancy is far enough from the dislocation, the interaction 
energy between both defects can be predicted using the elasticity theory 
Considering only the first-order size interaction, this energy is simply given 
by [1,2,29,30] 

= P sny, (2) 

where P = — (J^u <^ii) /3 is the pressure created by the dislocations at the lat- 
tice site where the vacancy lies and SQy is the vacancy relaxation volume. This 
last term is calculated directly from the atomic potential whereas the pressure 
is computed within the framework of the elasticity theory of dislocations [30]. 
Usually such calculations are performed assuming perfect dislocations and 
isotropic elasticity. However, dislocations are dissociated in fee metals and not 
all materials can be assumed isotropic. We therefore check the validity of these 
different simplifying assumptions. For the isotropic elasticity calculations, we 
use the shear modulus /x, the bulk modulus K and the Poisson coefficient u, 
obtained by Voigt average [30] of the elastic constants corresponding to the 
potential used for atomic calculations: 

1 (Cii - Ci2 + 3Cu) , (3a) 
o 

^ (Cn + 2C^2) , (3b) 
Cii + 4Ci2 — 2C44 ^^^^ 

2 (2Cii -|- 3Ci2 -l- C44) 



K = 



Whatever the assumptions used to calculate the pressure within the elastic- 
ity theory of dislocations, one has to cope with the problem of the periodic 
boundary conditions corresponding to the atomic simulations. Indeed, the two 
dislocations in the simulation box are not the only ones to consider in the elas- 
tic calculation. All their periodic images forming an infinite array of dipoles 
have to be included too. As the stress created by a dislocation is varying as the 
inverse of the distance, the sum associated with the calculation of the pres- 
sure is only conditionally convergent. To circumvent this numerical problem, 
we use the computational approach proposed by Cai^ [33,34]. 

^ So as to compute the stress field created by this double periodic array of dis- 
locations, one can take advantage too from the fact that each dislocation and its 
periodic images in the z direction form a tilt boundary. This causes a stress field 
which expression is analytical [30] and is known to decay exponentially with the 
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3.1 Perfect dislocations within the isotropic elasticity theory 



Assuming isotropic elasticity, the pressure created by a perfect edge dislocation 
at a point of coordinates {x,y,z)'^ is 



St: 1 — u + z'^ 



where Si = +1 (conversely —1) if the two supplementary half- plane corre- 
sponding to the edge dislocation i is in the z > (conversely 2; < 0) sub-space 
and b — a\/2/2 is the Burgers vector of the perfect edge dislocation. This 
pressure, predicted by isotropic elasticity, can be used in Eq. 2 after summing 
pressures arising from all image dislocations so as to predict the interaction 
energy between the vacancy and the dislocations. In Fig. 5, we compare the 
elasticity theory predictions with values of Ey^^l, previously calculated with 
EAM potentials, for different separation distances between both defects in the 
[110] and [111] directions. The elasticity theory manages to follow the atomic 
calculations as long as the vacancy is not too close to the dislocation core. For 
Al, a reasonable agreement is obtained for distances greater than ~ 10 A in 
[111] direction and greater than ~ 20 A in [110] direction. We believe that the 
elasticity theory is less accurate when the pressure is close to zero. 



distance normal to the wall. 

^ The origin is defined as the center of the perfect dislocation and axes are defined 
on Fig. la. 




Distance (A) Distance (A) 



(a) [110] direction. (b) [111] direction. 

Fig. 5. Variation with the separation distance of the vacancy - dislocation interaction 
energy Ey^^ in [110] and [111] directions in Al. Symbols correspond to atomic 
simulation and lines to the isotropic elasticity theory (Eq. 2) assuming perfect or 
dissociated dislocations. 
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3.2 Influence of dislocations dissociation 



We then examine whether considering the dissociation of the dislocations 
modifies the value of the interaction energy with the vacancy. According to 
isotropic elasticity theory, the dissociation distance is [29, 30] 

" 247771-1.' 

As presented in table 1, the dissociation distance predicted is quite close to 
the one observed in atomic simulations. 

Only the edge components of each partial dislocations contribute to the pres- 
sure which bewomes 

fib 1 + u / z z \ , , 

' " '6^1^ \{x-d^-y + z^ ^ {x + d'^^y + z^ J' ^ ' 

When the vacancy is close to the dislocation, the interaction energy £'v-d 
calculated using the pressure, as given by Eq. 6, is in better agreement with 
atomic calculations than with the pressure given by Eq. 4. Even for Al, which 
has the smallest dissociation distance, the improvement is obvious (Fig. 5). 
When the separation distance between the vacancy and the dislocation be- 
comes large compared to the dissociation distance {^/x^ + z^ » the 
difference between the dissociated and the non-dissociated cases becomes neg- 
ligible since Eq. 6 converges to Eq. 4. 



3.3 Influence of anisotropy 



Instead of assuming isotropic elasticity to calculate the pressure created by a 
dislocation, one can use the elastic constants Cn, C12 and C44 corresponding 
to the empirical potential of the atomic simulations and take full account of 
the anisotropy. Considering anisotropy can be quite significant in elastic cal- 
culations. For instance, two dilatation centers do not interact within isotropic 
elasticity, whereas there is an interaction when anisotropy is introduced. More- 
over, Cai [33,34] showed that anisotropy has to be considered when calculating 
dislocation basic properties like their core energy or the Peierls stress. 

Due to the symmetry in the [112] direction along the dislocation line, finding 
the stress field associated with an edge dislocation is a two-dimensional prob- 
lem and can be solved numerically following the sextic anisotropic elasticity 
theory of straight dislocations [30]. We thus calculate the elastic interaction 
between the two partial dislocations composing an edge dislocation so as to 
obtain its dissociation distance The obtained value is really close to the 
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one calculated with isotropic elasticity (Tab. 1). The agreement with the dis- 
tance observed in atomic simulation is slightly better but the difference is not 
really relevant, even for Cu which is the most anisotropic metal considered. 

Using anisotropic elasticity theory to compute the pressure created by the 
dissociated edge dislocations, Eq. 2 gives another estimate for the interaction 
energy with vacancies. As shown for Cu in Fig. 6, the values of Ey^l, are 
similar to the ones calculated with isotropic elasticity. The agreement with 
atomic calculations is slightly better, but, as for the dissociation distance, 
the change due to anisotropy is negligible. Anisotropy does not lead to any 
significant change compared to isotropic calculations. Therefore, it does not 
need to be considered in the elastic calculation of the interaction energy of a 
vacancy with dislocations. 



3.4 Inhomogeneity interaction 



Another elastic interaction between the vacancy and the dislocation, the in- 
homogeneity contribution, needs to be considered [1,5]. It arises from the fact 
that the vacancy has different elastic constants from the matrix. Thus, there 
is a change in the elastic energy that can be stored in the material. As the 
elastic constants that can be associated to the vacancy have to be lower than 
the ones of the matrix, this inhomogeneity interaction is always attractive. 
Indeed, the presence of a vacancy releases some elastic energy. Assuming that 
the elastic constants of the vacancy are null, one maximizes this contribution. 
In isotropic elasticity, the interaction energy between the vacancy and the 




-100 -50 50 100 -60 -40 -20 20 40 60 

Distance (A) Distance (A) 



(a) [110] direction. (b) [111] direction. 

Fig. 6. Variation with the separation distance of the vacancy - dislocation interaction 
energy ii^v-D [HO] [IH] directions in Cu. Symbols correspond to atomic 
simulation and lines to the isotropic or anisotropic elasticity theory (Eq. 2) assuming 
dissociated dislocations. 



13 



dislocations now reads 



Pinter _ p ^ Q ( ^(l - ^) 2 , 5(1 - u) 2] (^^ 

where a is the Von-Mises equivalent stress and is related to the pure shear 
components of the stress created by the dislocation by 

^EK + ^M'- (8) 

Like the pressure, it can be easily computed for dissociated edge dislocations 
using isotropic elasticity. Considering the inhomogeneity contribution, the in- 
teraction energy now depends not only on the pressure but also on the shear 
stress. As it is varying linearly with the square of the stress, and thus with the 
square of the inverse of the distance between both defects, this contribution 
is much shorter range than the size contribution. 

Considering the inhomogeneity contribution in the isotropic elasticity calcu- 
lation of EY^fji improves the agreement with atomic simulations close to the 
core of the partial dislocations, as can be seen for Ni in Fig. 7. In particu- 
lar, without this contribution the interaction is symmetric with respect to the 
dislocation glide plane. The inhomogeneity interaction thus explains why the 
vacancy attraction in the compressive region is larger in magnitude than its 
repulsion in the tensile zone. This non-symmetric behavior of the interaction 
energy is clearly visible when the vacancy lies no more than one atomic dis- 
tance from the dislocation core (Fig. 7(b)). For larger separation distances, 
the change induced by the inhomogeneity contribution is not significant. As 




-60 -40 -20 20 40 60 ' -60 ^0 -20 20 40 60 

Distance (A) Distance (A) 



(a) [110] direction. (b) [111] direction. 

Fig. 7. Variation with the separation distance of the vacancy - dislocation interac- 
tion energy E^^^^ in [110] and [111] directions in Ni. Symbols correspond to atomic 
simulation and lines to the isotropic elasticity theory assuming dissociated disloca- 
tions. The first order approximation corresponds to Eq. 2 and the second order to 
Eq. 7. 
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pointed by BuUough and Newman [1], when this inhomogeneity interaction 
is taken into account in the elastic calculation, one should incorporate also 
the second-order size interaction arising from the non-linear elastic properties 
of the crystal, as both contributions are decreasing in the same manner with 
the distance. Eventhough this last term may lead to a better agreement with 
atomic calculations in the vicinity of the dislocation cores, we do not consider 
it since results based on non-linear elastic calculations are too complicated to 
handle. Therefore, we prefer to use only the linear elasticity theory, consider- 
ing the fact that one can switch to atomic calculations for short separation 
distances where the theory breaks down. 



4 Conclusions 

The interaction between vacancies and edge dislocations has been modeled in 
fee metals at different scales. When only a few atomic distances he between 
the vacancy and the cores of the partial dislocations or the stacking fault an 
atomic description of this interaction is needed. Using empirical potentials, 
one finds for fee metals that the most attractive position of the vacancy lies 
at the edges of the supplementary half-planes corresponding to the partial 
dislocations. Positions above the slip planes, in the compressive region, are 
always attractive whereas, below the slip planes, positions are repulsive or, 
for some metals like Cu and Ni, slightly attractive in the habit plane of the 
stacking fault. 

When the separation distance between both defects is greater, the elasticity 
theory manages to predict quantitatively the interaction energy. One does 
not need to take into account anisotropy in the calculations as this does not 
lead to any significant change compared to an isotropic calculation. On the 
other hand, better agreement with atomic simulations is obtained for short 
separation distances, when the dislocation dissociation is considered. Taking 
into account the inhomogeneity interaction and not only the size effects im- 
proves predictions of the elastic calculation too. The most significant change 
is that the interaction energy is not anymore symmetric with respect to the 
dislocation glide plane. 
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